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Abstract 

The scenario that some first-order phase transitions may have taken place in the early 
Universe offers us one of the most intriguing and fascinating questions in cosmology. 
Indeed, the role played by the latent "heat" or energy released in the phase transition is 
highly nontrivial and may lead to some surprising, important results. In this paper, we 
take the wisdom that the cosmological QCD phase transition, which happened at a time 
between 10~^ sec and 10"'' sec or at the temperature of about 150 MeV and accounts 
for confinement of quarks and gluons to within hadrons, would be of first order. To get 
the essence out of the scenario, it is sufficient to approximate the true QCD vacuum 
as one of degenerate 0-vacua and when necessary we try to model it effectively via a 
complex scalar field with spontaneous symmetry breaking. We examine how and when 
"pasted" or "patched" domain walls are formed, how long such walls evolve in the long 
run, and we believe that the significant portion of dark matter could be accounted for in 
terms of such domain- wall structure and its remnants. Of course, the cosmological QCD 
phase transition happened in the way such that the false vacua associated with baryons 
and many other color-singlet objects did not disappear (that is, using the bag-model 
language, there are bags of radius 1.0 fcrmi for the baryons) - but the amount of the 
energy remained in the false vacua is negligible. The latent energy released due to the 
conversion of the false vacua to the true vacua, in the form of " pasted" or " patched" 
domain walls in the short run and their numerous evolved objects, should make the 
concept of the "radiation-dominated" epoch, or of the "matter-dominated" epoch to be 
re-examined. 

PACS Indices: 12.39.Ba, 12.38.Mh, 12.38.Lg, 98.80.Bp, 98.80.-k. 

1 Introduction 

The discovery |lj of fluctuations or anisotropies, at the level of 10~^, associated with the 
cosmic microvi^ave background (CMB) has helped transformed the physics of the early uni- 
verse into a main-stream research area in astronomy and in particle astrophysics, both 
theoretically and observationally [2j. CMB anisotropies [3] and polarizations [1], the latter 
even smaller and at the level of 10""^, either primary (as imprinted on the last scattering 
surface just before the universe was (379 it 8) x 10^ years old) or secondary (as might be 
caused by the interactions of CMB photons with large-scale structures along the line of 
sight), are linked closely to the inhomogeneities produced in the early universe. 

Over the last three decades, the standard model of particle physics has been well estab- 
lished to the precision level of 10~^ or better in the electroweak sector, or to the level of 
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10~'^ — 10~^ for the strong interactions. In the theory, the electroweak (EW) phase transi- 
tion, which endows masses to the various particles, and the QCD phase transition, which 
gives rise to confinement of quarks and gluons within hadrons in the true QCD vacuum, 
are two well-estabhshed phenomena. Presumably, the EW and QCD phase transitions 
would have taken place in the early universe, respectively, at around 10^^^ sec and at a 
time between 10~^ sec and 10~^ sec, or at the temperature of about 300 GeV and of about 
150 Mey, respectively. Indeed, it has become imperative to formulate the EW and QCD 
phase transitions in the early universe if a quantitative theory of cosmology can ever be 
reached. 

The purpose of this work is to focus our attention on cosmological QCD phase tran- 
sition and to assess whether its roles in the early universe can be synthesized in a more 
quantitative terms - a problem which has become one of the most challenging problems in 
the physics of the early universe. To simplify the situation, we use the bag-model language 
and try to model the degenerate 0-vacua, the lower-temperature phase, as the minima of 
the spontaneously-broken complex scalar fields. In particular, we try to set up the strategy 
of tackling the problem by dividing it into problems in four different categories, viz.: (1) 
how a bubble of different vacuum grows or shrinks; (2) how two growing bubbles collide 
or squeeze, and merging, with each other; (3) how the Universe eventually stabilize itself 
later while keeping expanding for several orders of magnitude; and (4) how specific objects, 
such as back holes or magnetic strings, get produced during the specific phase transition. 
Questions related to part (4), which are quite complicated, will not be addressed here; see, 
e.g., ref. [5]. In the framework which we consider, we could describe the intermediate solu- 
tions based on the so-called "pasted" or "patched" domain walls when the majority of the 
false vacua get first eliminated - but how it would evolve from there and how long it would 
evolve still uncertain. 

The major result of this paper is that the latent heat (or latent energy), 
which turns out to be identified as the "bag constant", is huge compared to 
the radiation density at the cosmological QCD phase transition (i.e. at about 
3 X lO^^sec). As time evolved to the present, the percentage of this quantity 
becomes probably the majority of dark matter (25 % of the composition of the 
present Universe). 

2 The Background Universe as from Einstein's General Rel- 
ativity and the Cosmological Principle 

A prevailing view regarding our universe is that it originates from the joint making of 
Einstein's general relativity and the cosmological principle while the observed anisotropics 
associated with the cosmic microwave background (CMB), at the level of about one part in 
100,000, may stem, e.g., from quantum fluctuations in the inflation era. In what follows, 
we wish to first outline very briefly a few key points in the standard scenario so that we 
shall have a framework which we may employ to elucidate the roles of phase transitions in 
the early universe. 

Based upon the cosmological principle which state that our universe is homogeneous 
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and isotropic, we use the Robertson- Walker metric to describe our universe[6]. 

ds^ = - R\t){^^ + r^dO^ + r^sin'ed^''}. (1) 
1 — Kr^ 

Here the parameter k describes the spatial curvature with k = +1, —1, and referring to 
an open, closed, and flat universe, respectively. The scale factor R{t) describes the size of 
the universe at time t. 

To a reasonable first approximation, the universe can be described by a perfect fluid, 
i.e., a fluid with the energy-momentum tensor T^^ y = diag {p, , —p, —p, —p) where p is the 
energy density and p the pressure. Thus, the Einstein equation, y = SttGnT^ y + Ag^ y, 
gives rise to only two independent equations, i.e., from {fi, u) = (0, 0) and (i, i) components, 

i?2 k SttGn a 

Combining with the equation of state (EOS), i.e. the relation between the pressure p and 
the energy density p, we can solve the three functions i2(t), p{t), and p{t) from the three 
equations. Further, the above two equations yields 

R A-kGn , ^ A 

showing either that there is a positive cosmological constant or that p+'ip must be somehow 
negative, if the major conclusion of the Supernovae Cosmology Project is correct |7j, i.e. 
the expansion of our universe still accelerating [R/R > 0). 

Assuming a simple equation of state, p = wp, we obtain, from Eqs. (2) and (3), 

2^ + (1 + 3w){^ + ^) _ (1 + nj)A = 0, (5) 

which is applicable when a particular component dominates over the others - such as in the 
inflation era (before the hot big bang era), the radiation-dominated universe (e.g. the early 
stage of the hot big bang era), and the matter-dominated universe (i.e., the late stage of 
the hot big bang era, before the dark energy sets in to dominate everything else). In light 
of cosmological QCD phase transition, we would like to examine if the radiation- dominate 
universe and the matter- dominated universe could ever exist at all, since this has become a 
dogma in the thinking of our Universe. 

For the Inflation Era, we could write p = —p and /c = (for simplicity), so that 

i?-^ = 0, (6) 

which has an exponentially growing, or decaying, solution R oc e^"*, compatible with the 
so-called "inflation" or "big inflation". In fact, considering the simplest case of a real scalar 
field </>(t), we have 

P=\^'' + V{(^), p = ^4>^-V{(^), (7) 
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so that, when the "kinetic" term ^4>'^ is neghgible, we have an equation of state, p ~ —p. 
In addition to its possible role as the "inflaton" responsible for inflation, such field has also 
been invoked to explain the accelerating expansion of the present universe, as dubbed as 
" quintessence" or " complex quintessence" [8] . 

Let's look at the standard textbook argument leading to the radiation-dominated uni- 
verse and the matter-dominated universe: 

For the Radiation-Dominated Universe, we have p = p/3. For simplicity, we assume 
that the curvature is zero {k = 0) and that the cosmological constant is negligible (A = 0). 
In this case, we find from Eq. (5) 

i?cxti (8) 

Another simple consequence of the homogeneous model is to derive the continuity equa- 
tion from Eqs. (2) and (3): 

d{pR^) + pd{R^) = 0. (9) 

Accordingly, we have p oc for a radiation-dominated universe {p = p/3) while p oc R^^ 
for a matter-dominated universe (p << p). The present universe is believed to have a 
matter content of about 5%, or of the density of about 5 x lO^^^g/cm^, much bigger than its 
radiation content 5 x 10^^^(7/cm^, as estimated from the 3° black-body radiation. However, 
as i — > 0, we anticipate i? — > 0, extrapolated back to a very small universe as compared to 
the present one. Therefore, the universe is necessarily dominated by the radiation during 
its early enough epochs. 

For the radiation-dominated early epochs of the universe with k = and A = 0, we 
could deduce, also from Eqs. (2) and (3), 

T = { ^lQ^h-^/\^K). (10) 
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These equations tell us a few important times in the early universe, such as 10~^^sec when 
the temperature T is around 300 GeV during which the electroweak (EW) phase transition is 
expected to occur, or somewhere between 10~ sec 300 MeV) and lO-^ec 100 MeV) 
during which quarks and gluons undergo the QCD confinement phase transition. 

For the Matter-Dominated Universe, we have p ~ 0, together with the assumption 
that /c = and A = 0. Eq. (5) yields 

R(xtt (11) 

As mentioned earlier, the matter density pm scales like R~^, or pm oc the latter similar 
in the radiation-dominated case. 

When t = lO^sec, we have p^ = 6.4 x 10^^^ gm / cm^ and pm = 3.2 x lO^^^gm/cm^ , 
which are close to each other and it is almost near the end of the radiation-dominated 
universe. The present age of the Universe is 13.7 billion years - for a large part of it, it is 
matter-dominated although now we have plenty of dark energy (65% ~ 70%). 

However, it is generally believed that our present universe is already dominated by the 
dark energy (the simplest form being of the cosmological constant; about 70%) and the 
dark matter (about 25%). The question is when this was so - when the dark part became 
dominant. 
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There is another basic point - trivial but important. For both the electroweak and 
QCD phase transitions in the early Universe, if the phase transitions are described (approx- 
imately) by the complex fields ^, then the density of the system is given by 

p = + + Pm + (12) 

before or after or during the phase transition is being taking place. For the electroweak 
or QCD phase transition, we know that pm « P7, but the role played by is clearly of 
importance in our considerations. 

What would be missing in the standard textbook argument? We would come 

back in Section 6 to this important point, after we set up the general framework and have 
gained enough of insights. The crucial point is whether cosmological QCD phase transition 
is the first-order phase transition - if it is, there is the latent "heat" or energy released in the 
transition; the story would change dramatically if the amount of energy density turns out 
to he greater than either or p„i in the previous radiation- dominated or matter- dominated 
era. We have to question what happened if p^ would no longer be dominant in the radiation- 
dominated universe - the " new" dominant sector of the universe may not influence the " old" 
radiation-dominated piece but Einstein equation in principle no longer guarantee its validity 
(about the old, no- longer-dominant sector). 

3 The Cosmological QCD Phase Transition - the Big Picture 

Let's try to focus on the QCD phase transition in the early Universe, or on the cosmological 
QCD phase transition. 

At the temperature T > Tc ~ 150MeV, i.e., before the phase transition takes place, 
free quarks and gluons can roam anywhere. As the Universe expands and cools, eventually 
passing the critical temperature the bubbles nucleate here and there. These bubbles 
"explode", as we call it "exploding solitons". When it reaches the "supercooling" temper- 
ature, Tg, or something similar, the previous bubbles become too many and in fact most of 
them become touched each other - now the false vacua or "bubbles" of different kind (where 
quarks and gluons can move freely) start to collapse - or we call it "imploding solitons". 
When all these bubbles of different kind implode completely, the phase transition is now 
complete. 

There is some specialty regarding the QCD phase transition in the early Universe. 
Namely, the collapse of the false vacuum does depend on the inside quark-gluon content 
- e.g., if we have a three-quark color-singlet combination inside, the collapse of the false 
vacuum would stop (or stabilize) at a certain radius (we called the bag radius, like in the 
MIT bag radius); of course, there are meson configurations, glueballs, hybrids, six-quark or 
multi-quark configurations, etc. The QCD phase transition in the early Universe does not 
eliminate all the false vacua; rather, the end state of the transition could have at least lots 
of baryon or meson states, each of them has some false vacuum to stabilize the system. 

How big can a bubble grow? It is with the fastest speed which the bubble can grow is 
through the speed of light or close to the speed of light. The bubble could sustain from the 
moment it creates, say, T ~ Tc to the moment of supercooling, Tg 0.95 • T^, or during the 
time span f ~ 3 x 10~^ x 0.05sec (or 1.5 x 10~^sec). So, the bubble can at most grow into 
c • 1.5 X 10~^sec or 4.5 x 10'^ cm. 
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How big was the Universe during the cosmological QCD phase transition? Compared to 
the size now, an estimate is the expansion of 5.7 x 10^^, a huge factor. (See the beginning 
of Sect. 6.) In the simplest approximation (when the scalar fields don't couple to the other 
fields, such as gluons or quarks), the domain walls cannot disappear - not only sometime 
because of the possible nontrivial topology but that there should be some QCD dynamics 
to annihilate the walls. In light of the huge expansion factor, the domain wall structure 
cannot survive, except the strict topology (which we call "domain- wall nuggets"). 

As a yardstick, we note that, at f ~ 10~^ sec or T ~ 300 MeV, we have 

~ 6.4 X 10^°c/m/cm^ ~ 3.2 x lO^^m/cm^. (13) 

Or, at t ~ 3.30 x lO^^ sec or T = ~ IBOMeV, we have 

Py = 5.88 X lO^gm/cm^, pm = 6.51 x lO'^gm/cm^. (14) 

Slightly later when QCD phase transition has completed, at t ~ 10~^ sec or T 100 MeV, 
we have 

~ 6.4 X lO^gm/cm^, pm ~ 1.0 x lO'^gm/crr?. (15) 

In what follows, we use the so-called "bag models" [9l [10] to have the simplified version of 
quark confinement - I think it is important to use the simplified version in the complicated 
cosmological environment, in order to extract meaningful results. 

When the low-temperature bubbles start to show up (i.e. to nucleate), it is about 
Tc ^ T < Tg. This period is to be called "exploding soliton era",[IO] to be described in 
the next section (Sect. 4). The supercooling temperature Ts, presumably ~ 0.95rc (to be 
determined more reliably in the future), refers to the situation where the bubbles begin 
to saturate. We call it the "colliding soliton era". This is to be described in Section 5. 
Presumably it would be over when T < — (Tc — Tg) or longer. So, the cosmological QCD 
phase transition would be over when the Universe was 10~''sec old. The important things is 
that, because the phase transition is of the first order, it releases a huge amount of energy: 

pyac = 1-0163 X 10^'^gm/cm^, (16) 

clearly much bigger than the radiation density (cf. Eq. (14)). This quantity is in fact the 
same as "the zero-point energy". That is why we question the radiation-dominated 
universe. 

When the low-temperature bubbles fill up the space, the neighboring two bubbles would 
in general be labelled by different 9ij representing different but degenerate vacua - we as- 
sume that there are infinite many choices of 6; they are degenerate but complete equivalent. 
The domain wall is used to separate the two regions. Three different regions would meet 
in a line - which we call a vortex. We have to estimate the total energy associated with 
the domain walls and the vortices - particularly when these objects persist to live on for a 
"long" time - say, r ^ 10~^sec. These domain walls and vortices are governed, in the QCD 
phase transition in the early Universe, by the QCD dynamics - this is an important point; 
if not, what else? It is a tough question to figure out how long the Universe would stabilize 
itself through QCD dynamics and others; my rough guess is from a few seconds to years, 
say r. 
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For the moment, QCD enables us to make some estimates. Let us focus on f ~ 10~^sec, 
where pm = 1-0 x IQ"^ gm/cm?. Or, considering a unit volume of 1.0 cm^, the amount of 
the matter would be lOOym or 5.609 x IQ'^^GeV / . One proton or neutron weighs about 
1 GeV/(? so, in a volume 1.0 cm^ at t ~ 10~^sec, we had at least 5.609 x 10^-*^ baryons or, in 
the MIT bag model language, 5.609 x 10^^ bags ot R= 1.0 fermi false vacua associated with 
the system. To begin with, all the excited baryons and mesons, including topological objects, 
and multi-baryons, all have equal opportunities. But, remembering Icm^ = 10^^ fermi^, 
most space had to collapse into the true vacua with different 9i. 

4 Exploding Solitons 

We begin our study by examining an isolated bubble - expanding, that is, inside the bubble, 
it is the true vacuum labelled by some 6; outside the bubble, the false vacuum; we are 
thinking of the Universe cooling down and expand. That is, how the bubble nucleates in 
the false vacuum (high-temperature). Remember that this happened in the period of time 
when Tc^T < Tg. 

Consider a spherical wall of radius R and thickness A separating the true vacuum inside 
from the false vacuum outside. The energy density difference of the vacua is S, the bag 
constant in the most simplified situation, and the energy r per unit area associated with 
the surface tension on the separating wall is a quantity to be calculated but nevertheless is 
small compared to the latent heat. If the wall expands outward for a distance 5R, then the 
energy budget arising from the vacuum change is 

B ■ 47ri?2 -SR-T- A7r{{R + 6Rf - R^} = -p6V, (17) 

where p is the pressure and is so defined that a negative pressure would push the wall 
outward. (We use the notation r here, since a and p are reserved for other purposes.) 

When the surface tension energy required for making the wall bigger is much less than 
the latent heat required from the expansion of the bubble, the bubble of the stable vacuum 
inside will grow in an accelerating way, possibly resulting in explosive growth of the bubble. 
The scenario may be as follows: When the universe expands and cools, to a temperature 
slightly above the critical temperature Tg, bubbles of lower vacua will nucleate at the spots 
where either the temperature is lower, and lower than Tc, or the density is higher, and 
higher than the critical density pc- As the universe continues to expand and cool further, 
most places in the universe have the temperature slightly below T^, that is, the destiny 
arising from eternal expansion of the universe is driving the average temperature of the 
entire universe toward below the critical temperature. The universe must find a way to 
convert itself entirely into another vacuum, the true vacuum at the lower temperature. 

Therefore, we have a situation in which bubbles of true vacua pop up (nucleate) here 
and there, now and then, and each of them may grow explosively in the environment made 
of the false vacuum for now, but previously the true vacuum when the temperature was still 
near the critical temperature Tc. In the expanding universe which cooled down relatively 
rapidly, i.e. from Tc to the supercooling temperature Tg, the situation is awfully complicated. 
When the temperature becomes lower than Tg, the problem can be modelled, in the simplest 
way, by characterizing the vacuum structure by a complex scalar field interacting via the 
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potential V{(f)): 

V{cf>) = ^cl>*cl>+^{cf>*cf>f, /x2<0, A>0. (18) 

For T > Tc, we have iJ^'^iT) > and A > 0, so it is between and Tg when the situations 
are awfully complicated (and we try to avoid in this paper). Note also that, in the complex 
scalar field description, the true vacua have degeneracy described by a continuous real 
parameter 9. (j) = everywhere in the spacetime describes the false vacuum for the universe 
at a temperature below the critical temperature T^.. Consider the solution for a bubble 
of true vacuum in this environment. It is required that the field must satisfy the field 
equation everywhere in spacetime, including crossing the wall of thickness A to connect 
smoothly the true vacuum inside and the false vacuum outside. This is why we may call the 
bubble solution "a soliton", in the sense of a nontopological soliton of T.D. Lee's. However, 
the soliton grows in an accelerating way, or the name "exploding soliton". 

The situation must have changed so explosively that at a very short instant later the 
universe expands even further and cools to even a little more farther away from T,. and most 
places in the universe must be in the true vacuum, making the previously false vacuum 
shrink and fractured into small regions of false vacua, presumably dominantly in spherical 
shape, which is shrinking in an accelerating way, or " implosively" . Using again the complex 
scalar field as our language, we then have "imploding solitons". 

In what follows, we attempt to solve the problem of an exploding soliton, assuming that 
the values of both the potential parameters and A are fairly stable during the period of 
the soliton expansion. The scalar field must satisfy: 

The radius of the soliton is R{t) while the thickness of the wall is A: 

A 

(f> = (pQ, for r < Ro + vt - —, 

= 0, for r>Ro + vt+^, (20) 

with R{t) = Ro + vt and v the radial expansion velocity of the soliton. 
We may write 

(l) = f{r + vt); w = {l-v'^)r, (21) 
so that the field equation becomes 

fi + -f = a-v')-'Xf{\f\'-cl>l). (22) 

aw"^ w aw 

We will be looking for a solution of / across the wall so that it connects smoothly the 
true-vacuum solution inside and the false vacuum solution outside. 
Introducing g = wf{w), we find 

g" = (1 - v')-'Xg{\ £ p -02}, (23) 

an equation which we may solve in exactly the same manner as the coUiding-wall problem 
to be elucidated in the next section. 
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5 Colliding Walls - Formation of " Pasted" Domain Walls 



When bubbles of true vacua grow explosively, the nearby pair of bubbles will soon squeeze 
or collide with each other, resulting in merging of the two bubbles while producing cosmo- 
logical objects that have specific coupling to the system. The situation is again extremely 
complicated. Remember that this happened when T ^ Tg, not too long after. 

We try to disentangle the complexities by looking at between the two bubble walls that 
are almost ready to touch and for the initial attempt neglecting the coupling of the vacuum 
dynamics to the matter content. Between the two bubble walls, especially between the 
centers of the two bubbles, it looks like a problem of plane walls in collision - and this is 
where we try to solve the problem to begin with. 

In fact, we have to consider one bubble first - the spherical situation as in the previous 
section but the bubble is "very" large we could look at the z-direction in the sufficiently 
good plane approximation (i.e. all bubble surfaces are just like planes). At this point, we 
have one wall, with thickness A, moving with velocity v in the z direction; on the left of 
the wall is the false vacuum, and on the right the true vacuum. 

The wall, of thickness A, separates the true vacuum on one side from the false vacuum 
on the other side of the wall. For the sake of simplicity, the wall is assumed parallel to the 
(xj/)— plane and are infinite in both the x and y directions. In addition, at some instant the 
wall is defined between z = and z = ^ with the instantaneous velocity +v. 

For z > R+ ^ and all x and y, the complex scalar field (j) assumes 0o, a value of the 
true vacuum (the ground state). On the other hand, for z < —R — y and all x and y, the 
complex scalar field assumes i;^ = 0, the false vacuum. As indicated earlier, the field (p 
must satisfy the field equation everywhere in spacetime: 

We may write the wall on the right hand side but moving toward the left with the 
velocity v: 

(f> = f(^z- vt), for z-vt>(}, t< R/v. (25) 

so that 

(1 - t;2)/" = A/(| / |2 C7 ^1 00 |> 0. (26) 

In fact, we are interested in the situation that the function in Eq. (20) is complex: 

/ = ue'^ (27) 

so that, with A = A/(l — v'^), 

u" -u{e'f = \u{u'^ -u"^), (28) 

2u'e' + ue" = 0. (29) 
Integrating the second equation, we find 

u^e' = K, (30) 
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with K an integration constant. The equation for u is thus given by 

K 



u" = ^ + \u{u^ -a^), (31) 



w 



provided that the 6 function is defined (in the region of the true vacuum and the wall). 

Let us try to focus on the last two basic equations - for u and 9, say, as the functions 
of ^ (e.g. ^ = z ±vt). For ^ > A, we have = ae^^ (the true vacuum) and, for ^ < 0, we 
have = (the false vacuum; with 9 undetermined). We find, for ^ ^ 0"*", 
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9 = - V^(/0(1 + F{0) + Co, (32) 

with Co a constant and i^(0 regular near ^ 0. Therefore the 9{^) function could be 
"mildly singular" or blow up near ^ ~ - this is in fact a very important point. 

Of course, the equation for u can be integrated out to obtain the result. For the "wall" 
region (i.e. < ^ < A), the solution reads as follows: 



0-2 ru^a^ fjy 

Jo ^-K + ay- 2/3y2 + j3y^ 



with 



(33) 

V — -fv -r ucy — ^/jy- -r jjy^ 
/•! dy 



^-K + ay- 2/?y2 + ^y3 



(34) 



Here (3 = ^cr^, and K and a parameters related to the integration constants. Of course, 
the solution in true- vacuum region can be obtained by extension. 

In the wall region, we could compute the surface energy per unit area (i.e. surface 
tension mentioned earlier in Eq. (16)): 



+ (35) 



some integral easy to calculate. 

There is an important note - that is, the solution for (f> obtained so far applies for the true 
vacuum and the wall, and which is continuous in the region; how about the false vacuum? 
This is an important question because in the false vacuum we know that u = but 9 is left 
undetermined. So, in first-order phase transitions we have certain function undefined in the 
false- vacuum region (s). This is a crucial point to keep in mind with. 

As a parenthetical footnote, we note that the equation for the exploding or imploding 
spherical soliton, Eq. (22), may be integrated and solved in an identical manner. 

Now let us focus on the merge of the two bubbles - the growing of the two true-vacuum 
bubbles such that the false-vacuum region gets squeezed away. This is another difficult 
dynamical question. In fact, we can make the false- vacuum region approaching to zero, 
i.e., the region with the solution u = gets squeezed away; one true- vacuum region with 
9i and Ai (the latter for the wall) is connected with the one with 92 and A2 - we could 
use (Ki,K2) to label the new boundary; to be precise, wc could call it "the pasted domain 
wall" or "the patched domain wall". It is in fact two walls pasted together - if we look at 
the boundary condition in between, we realize that the structure would persist there for 
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a while to go. The pasted domain wall could evolve further but this may not be relevant 
for counting the energies involved. The evolved forms of the pasted domain walls could be 
determined by the topology involved - for the purpose of this paper, we can ignore this fine 
aspect. 

Suppose that the cosmological QCD phase transition was just completed - we have 
to caution that, not everywhere, the false vacua be replaced by the true vacua so that in 
between the walls be replaced (approximately) by the pasted domain walls. There are places 
for color-singlet objects (i.e. hadrons) which quarks and gluons tried to hide; these places 
are still called by the "false vacua" with the volume energies. Thus, the volume energy, 
i.e. B in Eq. (16) or defined suitably via A and (in Eq. (17)), or at least some portion 
of it, may convert itself into the surface energy and others - B = 57 MeV/ fm^ using the 
so-called " bag constant" in the MIT bag model [9j or Columbia bag model [10] . 

This energy density B = blMeV/ fm? = 1.0163 x is huge as compared to 

the radiation density (which is much bigger than the matter density pm) at that time, 
t ~ 10~^ ~ 10~^sec (see Eqs. (13)- (15)). Some exercise indicates that this quantity of 
energy is exactly the latent "heat" or energy released in the first-order phase transition. 

The cosmological QCD phase transition should leave its QCD mark here - since the 
volume energy that stays with the "false vacuum" is simply reduced because the volumes 
with the "false vacua" are greatly reduced - but not eliminated because quarks and gluons, 
those objects with colors, still have some places to go (or, to hide themselves). 

6 Possible Connection with the Dark Matter 

Let us begin by making a simple estimate - the expansion factor since the QCD phase 
transition up to now. The present age of the Universe is 13.7 billion years or 13.7 x 10^ x 
365.25 X 24 X 3600 or 4.323 x 10^"^ seconds. As indicated earher (cf. the end of Sec. 2), 
about the first lO^sec period of the hot big bang is previously-believed radiation-dominated. 
Consider the length l.Ofermi at i ~ 10~^sec, it will be expanded by a factor of 10^ up 
to t ~ lO^sec (radiation-dominated) and expanded further by another factor of 5.7 x 10^ 
until the present time - so, a total expansion factor of 5.7 x 10^^; changing a length of 
2 fermi at t ~ 10~^sec into a distance of 1 cm now. A proton presumably of i? = 1 fermi 
at i ~ 10~^sec should be more or less of the same size now; or, the bag constant or the 
energy associated with the false vacuum should remain the same. 

What would happen to the pasted or patched domain walls as formed during the cosmo- 
logical QCD phase transition? According to Eqs. (30) and (31) together with Eq. (32), we 
realize that the solutions in previously two different true- vacuum regions cannot be matched 
naturally - unless the K values match accidently. On the other hand, it is certain that the 
system cannot be stretched or over-stretched by such enormous factor, 10^2 or igi^. 

As we said earlier, at some point after the supercooling temperature T^, say, at Ts — 
X{Tc — Tg) (with A an unknown factor, presumably A S> 1), the system (the Universe) was 
temporarily stabilized since most of the pasted or patched domain walls had no where to 
go. Remember that all these happened in a matter of a fraction of 10~^,sec, as judging from 
the size of Tc and Tg. The next thing to happen is probably the following. 

We believe that the field (p, being effective, cannot be lonely; that is, there are higher- 
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order interactions such as 

co(/'G;iG'^'^ ci(t>GGG, doHi', (36) 

some maybe being absent because of the nature of (p. In other words, we may beheve that 
the strong interactions are primarily responsible for the phase transition in question, such 
that the effective field (p couples to the gluon and quark fields; the details of the coupling 
are subject to investigations. 

That is, when the field (j) responsible for the pasted or patched domain walls is effective - 
the 4> field couples, in the higher-order (and thus weaker) sense, to the gluon and quark fields. 
It is very difficult to estimate what time is needed for pasted domain walls to disappear, if 
there are no nontrivial topology involved. If there is some sort of nontrivial topology present, 
there should left some kind of topological domain nugget - however, energy conservation 
should tell us that it cannot be expanded by too many orders (but our Universe did expand 
for many many orders of magnitude). I would guess that it takes about from a fraction of 
a second to several years (from the strong interaction nature of the problem), but certainly 
before the last scattering surface (i.e. 3.79 x 10^ years). 

To summarize, the energy associated with the cosmological QCD phase transition, 
mainly the vacuum energy associated with the false vacuum, disappeared in several ways, 
viz.: (1) the bag energies associated with the baryons and all the other color-singlet objects, 
(2) the energies with all kinds of topological domain nuggets or other topological objects, 
and (3) the decay products from pasted or patched domain walls with trivial topology. 

Let us begin with the critical temperature T = Tc ~ 150 MeV or t ~ 3.30 x 10~^sec. 
At this moment, we have 

p^ac = 1.0163 X W^'^gm/cm^, = 5.88 x W^gm/cm^, p„ = 6.51 x W'^gm/cm^. (37) 

Here the first term is what we expect the system to release - the so-called "latent heat"; I 
call it "latent energy" for obvious reasons. The identification of the latent "heat" with the 
bag constant is well-known in Coulomb bag models |10j . 

This can be considered just before the cosmological QCD phase transition which took 
place - at the moment the energy components which we should take into consideration. 

As time went on, the Universe expanded and the temperature cooled further - from the 
critical temperature to the supercooling temperature (T^ ~ 0.95 x Tc with the fraction 0.95 
in fact unknown) and even lower, and then the cosmological QCD phase transition was 
complete. When the phase transition was complete, we should estimate how the energy 
Pyac is to be divided. 

Let's assume that the QCD phase transition was completed at the point Tg (in fact 
maybe a little short after Tg). Let's take Tg = 0.95 Tc for simplicity. We would like to 
know how the energy p^ac is to be divided. First, we can estimate those remained with 
the baryons and other color-singlet objects - the lower limit is given by the estimate on the 
baryon number density (noting that one baryon weighs about l.OGeV/c^): 

Pm = 6.51 X lO'^gm/cm^ x 0.5609 x lO^^GeV / / gm = 3.65 x IQ^^GeV/ / cm^ . (38) 

So, in the volume l.Ocm^ or lO'^^/ermi^, we have at least 3.65 x lO'^^ baryons. One baryon 
has the volume energy (i.e. the bag energy or the false vacuum energy) 57MeV/fermi^ x 
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|7r(1.0/ermi)^ (which is 238.8MeF). So, in the volume l.Ocm^, we have at least 238.8Mey x 
3.65 X 10^6 or 8.72x lO^^GeF in baryon bag energy. Or, in different units 8. 72 x 1025/(0.5609 x 
10^^) gm/c^ or Ibh.fygm/ cP' . Only a tiny fraction of p^ac is to he hidden in baryons or other 
color-singlet objects after the QCD phase transition in the early Universe. 

So, where did the huge amount of the energy pyac go? In the beginning of the end of 
the phase transition, the pasted domain walls with the huge kinetic energies seem to be the 
main story. A pasted domain wall is forming by colliding two domain walls while eliminating 
the false vacuum in between. The kinetic energies associated with the previously head-on 
collision become vibration, center-of-mass motion, etc. Of course, the pasted domain walls 
would evolve much further such as through the decaying interactions given earlier or forming 
the "permanent" structures. In any case, the total energy involved is known reasonably - a 
large fraction of pyac, much larger than the radiation pj (with p^ negligible at this point). 

The story is relatively simple when the cosmological QCD phase transition was just 
completed and most "pasted" domain walls still have no time to evolve. We return to Eqs. 
(2) and (3) (i.e. Einstein equations) for the master equations together with the equation of 
state with p and p determined by the energy-momentum tensor: 



dL 

Tt. = 9i^Q(Q^f-^ - (39) 



Further analysis indicates that the equation of state for the "pasted" or "patched" domain 
walls is nothing unusual - the reason is that we are working in the real four-dimensional 
space-time and all of the objects are of finite dimensions in all the directions. The "domain 
walls" discussed by us are for real and cannot be stretched to infinity in a certain dimension. 

In fact, there is certain rule which one cannot escape. Let assume a simple equation of 
state, p = wp, for simplicity and come to look at Eq. (5). Let's consider the situation in 
which there is no curvature k = and the cosmological constant A is not yet important. 

2| + (1 + 3«;)|^ = 0, (40) 

which yields 

i? oc (41) 

From the equation of continuity, d{pR) + pd{R^) = 0, it is easy to obtain p oc Ep^^~^^\ 
Thus, we deduce that, under very general situations, the density behaves like 

p = Ct-\ (42) 

where the constant C is related to w in the simplified equation of state. It is clear that the 
limit to w = —1 (the cosmological constant) is a discontinuity. 
Of course, Eq. (4) is still valid: 

R^ 3~^^ 3"' 

This has an important consequence - the idea of the previous universe expansion usually 
based on the radiation alone from t ~ 10"^'' sec (after the cosmological electroweak phase 
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transition had taken place) to t ^ 10^ sec (when it was close that p-y = pm) has to be 
modified because the latent energy p^ac was about 2 x 10^ times the radiation energy at the 
moment of the cosmological QCD phase transition. 

Shown in Fig. 1 is our main result - even though it is a qualitative figure but it tells 
us a lot. At t 3.30 x 10~^ sec, where did the latent energy lO^^gm/cm^ evolve into? We 
should know that the curve for p^, for massless relativistical particles, is the steepest in 
slope. The other curve for p^ is the other limit for matter (which P « 0). In this way, the 
latent energy is connected naturally with the curve for pdm - in fact, there seems to be no 
other choice. Remember that p oc except the slope for different types of " matter" . 



Figure 1: The various densities of our universe versus time. 

Coming back to Eq. (43) or (4) , we could assume for simplicity that when the cosmolog- 
ical QCD just took place the system follows with the relativistical pace (i.e. P = p/2>) but 
when the system over-stretched enough and had evolved long enough it was diluted enough 
and became non-relativistic (i.e. P 0). It so happens that in both cases the density to 
the governing equation, Eq. (43) or (4), looks like p oc t^^ although it is i? oc followed 

2 

by i? OC 1 3 . 

It is so accidental that what we call "the radiation-dominated universe" is in fact dom- 
inated by the latent energy from the cosmological QCD phase transition in the form of 
"pasted" or "patched" domain walls and the various evolved objects. In our case, the tran- 
sition into the "matter-dominated universe", which happened at a time slightly different 
from t ~ lO^sec, occurred when all the evolutions of the pasted domain walls ceased or 
stopped. In other words, it is NOT the transition into the "matter-dominated universe", 
as we used to think of. 

In fact, the way of thinking of the "dark matter", or the majority of it, turns out to be 
very natural. Otherwise, where did the 25% content of our universe come from? Of course, 
one could argue about the large amount of the cosmological QCD phase transition. We 
believe that the curves in Fig. 1 make a lot of sense. 
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Of course, one should ask what would happen before the cosmological QCD phase 
transition. It might not be the radiation-dominated. I believe that it opens up a lot of 
important and basic questions. 

7 Summary 

To sum up, we tried to illustrate how to describe the QCD phase transition in the early 
Universe, or the cosmological QCD phase transition. 

The scenario that some first-order phase transitions may have taken place in the early 
Universe offers us one of the most intriguing and fascinating questions in cosmology. In fact, 
the role played by the latent "heat" or energy released in the process is highly nontrivial. 

In this paper, I take the wisdom that the QCD phase transition, which happened at 
a time t ^ 3.30 x 10~^ sec or at the temperature of about 150 MeV and accounts for 
confinement of quarks and gluons to within hadrons in the true QCD vacuum, would be of 
first order. Thus, it is sufficient to approximate the true QCD vacuum as one of degenerate 
0-vacua and when necessary we try to model it effectively via a complex scalar field with 
spontaneous symmetry breaking. We examine how and how long "pasted" or "patched" 
domain walls were formed, how and how long such walls evolve further, and why the majority 
of dark matter might be accounted for in terms of these evolved objects. 

Our central result could be summarized by Fig. 1 together with the explanations. 
Mainly, we are afraid that the "radiation-dominated" epoch and the "matter-dominated" 
epoch, in the conventional sense, could not exist once the cosmological QCD phase transition 
took place. That also explains why there is the 25% dark-matter content, larger than the 
baryon content, in our present universe. 

Footnote: During the period which the paper is revised and refereed, some early version 
of this paper has been accepted for published in Modem Physics Letters A. 
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